function [position,alpha_ind ] = alphaindicator( pred,vpred,data,C )
%ALPHAINDICATIOR Summary of this function goes here
%   Detailed explanation goes here 

N = length(data);
alpha_ind = ones(1,N+1);
alpha_ind(1) = 0;
alpha_ind(2:N+1) = (pred(2:N+1)-data(1:N))./(vpred(2:N+1).^(1/2));
position = zeros(1,N);
i =2;
for i=2:N+1
    if (alpha_ind(i)>C)
        position(i) = 1;
    else
        position(i) = 0;
    end
end


end

